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SUMMARY 


Theoretically derived charts showing the profile-drag — thrust ratio 
are presented for a helicopter rotor operating in forward flight and 
having hinged rectangular blades with a linear twist of - 8 °. The charts 
showing the profile-drag characteristics of the rotor for various com- 
binations of pitch angle ^ ratio of thrust coefficient to solidity, and a 
parameter representing shaft power input, are presented for tip-speed 
ratios ranging from 0.05 to O.5O. In addition, the ratio of thrust coef 
ficient to solidity as a function of inflow ratio and blade pitch angle 
is presented in chart form. 

The charts of this paper differ from the rotor performance papers 
previously published by the National Advisory Committee for Aeronautics 
in that the theory on which the charts are based includes an approximate 
allowance for stall in the reversed-flow region and contains no small- 
angle assumptions regarding blade-section inflow angles and velocities. 
The charts of this paper are therefore considered more accurate than 
previous ones for fli^t conditions involving high inflow velocities 
and large regions of reversed velocity that may be encountered by high- 
performance helicopters. The assumption is made, however, that outside 
of the reversed- velocity region, the section angles of attack are small j 
thus the angles can be replaced by their sine. In addition, other than 
including an approximate allowance for stall in the reversed-velocity 
region, the charts do not Include stall and compressibility effects. 

The charts may be used to study the effects of design changes on 
rotor performance and to indicate optimvim performance conditions, as 
well as to estimate quickly rotor performance in forward flight. They 
are also usef\xl in obtaining inflow-ratio and pitch-angle values for use 
in calculating flapping coefficients and spanwlse loadings. The method 
of applying the charts to performance estimation is illustrated through 
sample calculation of a typical rotor-performance problem. 
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INTRODUCTION 


Equations vere presented in reference 1 from vhich the thrust^ the 
accelerating and decelerating torque^ and the profile-drag power of a 
hinged rotor operating at high tip-speed ratios and inflow angles could 
be calculated* Because the equations do not place any limitation on the 
magnitude of the inflow angle or on the rotor angle of attack^ they are 
considered more accurate than previous analyses when applied to high- 
speed helicopters and to certain types of convertible aircraft. This 
paper is an extension of reference 1 in that the equations of that refer- 
ence are used as the basis of a method for calculating the performance of 
lifting rotors over a wide range of operating conditions. 

Because the basic equations are lengthy^ the application of the 
method is considerably simplified by presenting the more lengthy equa- 
tions in the form of charts from which rotor performance can be quickly 
estimated. The charts cover operation at any rotor angle of attack at 
tip-speed ratios varying from 0.05 to O.5O for blades that are twisted 
-8^ (blade pitch angles at the tip 8^ lower than at the root). With the 
exception of an approximate allowance for stall in the re versed- velocity 
region^ the charts do not include stall and compressibility effects. 

Limit lines showing conditions for which blade angles of attack 
exceed specified values at given radial stations are included in the 
charts. These limit lines are useful in determining operating conditions 
at which stalling begins and for determining the limiting operating 
conditions . 


SYMBOLS 


slope of curve of section lift coefficient against section 
angle of attack^ per radian (assumed equal herein to 5-73) 

number of blades per rotor 


Cl 


rotor lift coefficient^ 


2 


rotor-shaft power coefficient. 


nR^p(fiR)5 


"^T 


rotor thrust coefficient^ 


T 





MCA TN 5325 


5 


c 


c 


e 
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blade section chords ft 


equivalent blade chord (weighted on thrust basis) ^ 


section profile-drag coefficient 



section lift coefficient 


helicopter parasite drag, lb 


rotor profile drag-lift ratio 


para site -drag area, 




, sq ft 


mass moment of inertia of blade about flapping hinge, slug-ft^ 


rotor lift , lb 

rotor-shaft power, ft-lb/sec 

shaft-power parameter, where P (in this ratio only) is equal 
to rotor-shaft power divided by velocity along flight path 
and is therefore also equal to drag force that could be 
overcome by shaft power at flight velocity 

blade radius measured from center of rotation, ft 

radial distance from center of rotation to blade element, ft 

rotor thrust , lb 

true airspeed of helicopter along flight path, fps 
induced velocity at rotor (always positive), fps 
helicopter gross weight, lb 

ratio of blade-element radius to rotor-blade radius, r/R 
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a rotor angle of attack| angle between axis of no feathering 

(that is^ axis about which there is no cyclic -pitch change) 
and plane perpendicular to flight path, positive when axis 
is inclined rearward, deg 

blade-element angle of attack, measured from line of zero lift, 
deg (when used in three-term drag polar in fig. l(b), is 

exp>ressed in radians) 

blade-element angle of attack at any radial position x and at 
^ ^ any blade azimuth angle ilf, deg; for example, 0)(270^) 

is blade-element angle of attack at tip of retreating blade 
at 270^ azimuth position 


4) (270^) blade-element angle of attack at radius at which 

' tangential velocity equals 0.4 tip speed and at 

270^ azimuth position, deg 

7 flight-path angle (positive in climb, negative in glide), deg 


0 blade-section pitch angle at 0.75 radius; angle between line 

*"^5 of zero lift of blade section and plane perpendicular to 

axis of no feathering, deg 

02_ difference between blade root and blade-tip pitch angles, posi- 

tive when tip angle is larger, deg 


A 


inflow ratio. 


V sin g - V 

m 


^ tip-speed ratio, 

nR 

p mass density of air, slugs /cu ft 

a rotor solidity, bc^/itR 

^ blade azimuth angle measured from downwind position in direc- 

tion of rotation, deg 

n rotor angular velocity, radians/sec 

Subscripts : 
c climb 


i 


induced 
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o profile 

p parasite 

V vertical component 


METHOD OF ANALYSIS 


The performance method presented herein utilizes the equations 
developed in reference 1 for blade-flapping coefficients^ rotor thrust^ 
torque, and profile-drag power and also, with some modifications, the 
energy performance analysis described in reference 2. Inasmuch as the 
performance method described herein is based on the equations developed 
in reference 1, the assumptions and limitations incorporated in the 
reference eqioations also apply to the perfonnance calculations. (The 
effects of the primary assumptions and limitations are discussed sub- 
sequently in the section entitled “Range of Application of Charts.”) 

In utilizing the equations of reference 1 to compute rotor perform- 
ance, mmierical substitutions regarding section lift and drag character- 
istics were made by using the values shown in figure 1. The values of 
Cl and c^^ below the stall given in the figure are representative of 

” semismooth” blades and are the same values used in the construction of 
the charts of reference 2. The values of c^ and c^^ above the stall 

are based to some extent on wind-timnel data presented in reference 3^ 
and are presented in figure 1 on the concept of a 360 ^ angle -of -attack 
range. This concept is useful in the analysis because the angle of attack 
in the reversed- velocity region can exceed l 80 ^. 

By following the procedure of reference 1, it was assumed that the 
thrust, torque, and power contributions of the reversed- velocity region 
could be approximated by using constant lift and drag coefficients cor- 
responding to a single representative section angle of attack. For each 
flight condition, the representative angle was computed at a radial sta- 
tion about one-third of the distance from the center of rotation to the 
outboard edge of the reversed-velocity region and at an aximuth angle of 
270 ^. The forces at this radial station were found to represent approxi- 
mately the average of the forces in the reversed-velocity region from 
plots of the radial distribution of the forces determined from step-by- 
step calculations for several sample cases. The values of ci and c^^ 

corresponding to the representative angle of attack were obtained from 
figure 1. Although some loncertainty as to the maximum value of c^^ in 

the 90 ^ angle-of -attack region exists, it was found that the use of a 
maximum value of 2.0, for example, instead of 1.6 had a negligible effect 
on the chart values over the range of applicability of the charts. 



6 


mCA TN 5325 


Fimdamental Performance Equation 

The power supplied at the rotor shaft of a helicopter is expended 
in overcoming the rotor profile-drag losses^ the induced drag losses, 
and the parasite-drag losses and in changing the potential energy of the 
aircraft in climb. The division of shaft power between the various 
sources can be written in coefficient form as 


Cp = + Cp^ + Cpp + Cp^ 


( 1 ) 


In presenting the relationship between Cp and Cp^ for various flight 

conditions in chart form, the resulting plots are greatly clarified if 
the power-coefficient ratios are divided by the thrust coefficient. Thus, 


Ct 


Ct 


+ 


Ct 


Ct Ct 


(2) 


Each ratio of power coefficient to thrust coefficient in equation ( 2 ) may 
be considered alternately as either an equivalent drag-thrust. ratio 
(wherein the equivalent drag is equal to the drag force that would absorb 
the power at a velocity equal to HR) or as an efficiency factor repre- 
senting power per unit thrust at a given tip speed. 

Almost any problem in helicopter performance, whether it be to deter- 
mine the shaft power required to maintain a steady-flight condition, the 
rate of climb at a given power condition, or the top speed of a given 
helicopter can be solved by means of the fundamental power relation 
expressed by equation ( 2 ). It will be noted that the familiar P/L, 
(d/l)q, . . . notation used in previous NACA helicopter performance papers 
is replaced herein by Cp/CQi, Cp^^C^, .... The power coefficients are 

based on the relatively constant HR instead of on V; thus, the nota- 
tion used herein avoids having the equivalent drag approach infinity as 
the forward speed approaches zero. For the same reason, rotor lift L — 
based on C^, which is dependent on forward speed — is replaced by the 

rotor thrust T inasmuch as Ctp is independent of forward speed. The 
conversion of one form of ratio to another is simply: 



(5) 



NACA TN 5525 


7 


Relations Required in Performance Uaxculations 

Formulas that are necessary for evaluating helicopter performance hy 
means of equation (2) are listed as follows; 

T cos(a +7) = W + Dp sin 7 (^) 



These equations^ with the exception of equation ( 9 )^ similar to those 
derived in chapter 9 of reference 4 except for the factor (i/cos a. Equa- 
tion (9) includes a drag term and is derived from a corresponding equation 
presented in the appendix of reference 5 "by using a miiltlplying factor of 
n/cos a. 
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Performance Charts 

The calculation of the various Cp/Ciji ratios in equations ( 5 ) to ( 9 ) 
can be greatly simplified by means of charts that relate the more lengthy 
ratios to the fimdamental variables 9 and |jl and to each other. 

Such charts are presented in figures 2 to and their use is demonstrated 
in succeeding sections of this paper. 


Each chart of figure 2 gives 2Gp/aa as a fiznction of A and Q .75 
for fixed values of \x ranging from 0.05 to O. 5 O. In figure 


is shown as a function of 


and 0 


.75 


for fixed values of 


^p/^T^ 2Cpycraj 

\i ranging from 0.05 to O. 5 O- Also, stall limit lines, the significance 
of which is discussed in references 2 and 6, are shown in these plots. 
Figure 4 is a graphical presentation of equation ( 9 ) from which the climb 
Cp 

parameter _S. may be determined from the climb angle 


n ’^T 
para site “drag parameter 


7 and the 


Cp. 


p cos a 


Gp 


OUTLINE OF PERPORMANCE METHOD UTILIZING CHARTS 


The problem of computing helicopter performance may be thought of 
as finding the value of one variable for given values of other pertinent 
variables, the variables being related by a number of basic equations. 
The problem, in essence, thus becomes the solution of a number of simul- 
taneous equations. The procedure can be greatly simplified by utilizing 
the performance charts presented in figures 2 to 4. The steps required 
in two typical types of performance calculations will be outlined and 
demonstrated by a sample calculation. 


Calculation of Rate-of -Climb Curves 

If the rate of climb (or descent) is required, the calculating pro- 
cedure would be as follows (for a given p) for the known parameters P, 
W, a, f, HR, and p: 

( 1 ) Assume T ~ W and calculate Cp. 

( 2 ) Compute Cp^Cp from equation (5)* 

( 3 ) Find Cp^^Cp and from figure 3 • 

(4) Find A from figure 2. 
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(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 


Calculate 

Compute 

Compute 


a from equation (lO) . 

Cpi^Cij amd Cp^JC^ from equations ( 7 ) and (8) . 
^Pc/ Crp f*poin ©(jusision (2) « 


Find 7 from figure 4. 

Compute Yy from the relationship = V sin y. 


(lO) If both 7 and Dp are very large, a new Crj\ can he computed 
by means of equation (4) and the process repeated to find a new value of 

Vv 


Calculation of Power-Required Ciirves 

A common performance calculation is to find the power required by 
a helicopter flying at a given airspeed and at a given rate of climb (or 
climb angle). The procedure would be as follows for the known param- 
eters W, a, f, 7 , p, and V: 

p 

(1) Assume that a = 0° and that (A/p) « 1; then, calculate T 

(and Qp) from equation (^). 

( 2 ) Calculate p from its definition. 

( 3 ) Calculate ^Pc/^ from equations ( 7 ), 

( 8 ), and ( 9 ), respectively. For convenience of application equation ( 9 ) 
has been used to construct figure from which can be obtained, for 
example, values of Cp^^C^ for given values of 7 and Cp^y^C^. 

(4) On the appropriate chart of figure 3^ lay off the sum of Cp ^C^, 

Cpp^Cp, and Cp^y^C^ along the Cp^Cip axis. Then, with that point as a 
base, construct a line having a slope of unity. (if ORq/^T ^p/^T 

were drawn to the same scale in figure 3? the construction line would be 
drawn at a 45^ angle. For the actual scales of figure 3^ the line is 
constructed at an angle which has a tangent of 2.) The intersection of 
this line with the proper line will yield values for 
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Cp^Cip, and 0 This procedure is illustrated for a constant tip-speed 

ratio by the following sketch: 



In order to avoid interpolation between tip-speed-ratio charts^ the 
value of V can be chosen so that ji is an even multiple of 0.05j 
otherwise^ the answer can be linearly interpolated between two successive 
charts . 


(5) Since 
from figure 2. 



0 and 


[1 are now known ^ 


X can be found 


(6) Compute a from equation (lO) and recompute \x from the equa- 


tion 


V cos g 

m 


(7) Recompute Cp^/Cijxy find new values 

of Cp/C^ and these values differ from the initially com- 

puted ones by more than a few percent^ repeat the process. Nonmlly one 
iteration is s-ufficient. However^ when a is within the range of ^20^ 
and M- ^ 0.50^ the initial assumptions that cos a = 1 and ("A/p) « 1 

are adequate and no iterations are needed. 
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Sample Performance Calciilation 


The performance calculations outlined in the preceding section will 
he illustrated by a sample problem; Calculate the power required by a 
helicopter traveling at I80 feet per second and climbing at a rate of 
500 feet per minute. The following additional data are known: 

W = 4,287 pounds, a = 0.08, fiR = 60O feet per second, p = 0.00238 slugs 
per cubic foot, R = 20 feet, 9 ^^ = -8°, and f = 12 square feet. 


(1) Assume that a = 0 ° and (A/p)^ « 1 . Also, 


7 = sin"-^ = sin "" = 1.6° and Dp = = 463 pounds. Then, from 


-1 - sin"^ 

T ' 180 ” “P - 2 

equation T = 4^500 pounds^ and Gp = O.OCA-0. 

( 2 ) Then, p = iQojSoo = 0.30. 

( 3 a) From equation ( 7 )^ ~ O.OOST. 

(31) From equation ( 8 ), Cp /C^ = 0 . 0322 . 

Cp ^ 

( 3 c) Then, — — = 0 . 107 - From figure 4 , 


P Cip 


cos a 
P 


Cp 


2 . = 0.028. 


"T 


Thus, ^ = 0 . 0084 . 


T 

Cp . Cp_ Cp 

(ka) + — 2. = 0.0067 + 0.0322 + 0.0084 = 0.0473. 

Cq> Ciji 

( 4 b) For p = 0.30 and 2Gji^cra = O.OI8, figure 3 (e) gives: 


CPr 


-/ip 


= 0.0315 


= 0.0788 


"T 


9 = 9 ^ 

.75 ^ 


(5) For 9 = 9 ° and 2 C^/aa = O.OI8, figure 2 (e) gives 


A = -0.080. 
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(6) The rotor angle of attack a can now be computed from equa- 
tion (10) as follows: 


^ -0.080 . 

tan a = + 

0.50 


0.00k 


0.18(1 + 0.0712) 


1/2 


-0.2k5 


a = -13.8° 

(7a) Recomputing the power coefficients with the above values for a 
and A results in changes that are within the accuracy of the computa- 
tions; therefore^ the originally computed values are sufficient. 

(7b) The power required is then calculated as 


p 

Power = — 

^T 

= ( 0 . 079) (0 . 00k )rt ( 20) ^ ( 0 . 00238 ) (600 ) ^ 

= 20k, 000 ft -lb /sec 

= 371 kp 

(7c) The rotor profile-drag power is 


Profile power = X 20k, 000 

0.079 

= 80,000 ft-lb/sec 
= lk6 hp 


RANGE OF APPLICATION OF CHARTS 


In the preparation of the charts, it was necessary to make some 
assumptions regarding the rotor physical parameters to be used with the 
theory. Some of the more pertinent effects of these assumptions as well 
as the effects of the restrictive assumptions of the theory are discussed 
in the succeeding sections. 
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Blade Characteristics 

The sample rotor for which the charts presented herein were prepared 
was assumed to have hinged rectangular blades with a mass factor pcaR^/lj_ 
equal to 15 and a linear twist of -8^. However^ according to the error 
analysis made in reference 6^ it would appear that the charts would be 
applicable to rotors having values of mass factors ranging from 0 to 25* 
Thus^ although blade-flapping motion is sensitive to mass factor, average 
rotor forces are relatively insensitive to moderate changes in the flapping 
motion. 

Although the charts were calculated for rotors having xmiform- chord 
blades, previous experience has shown that, in general, the forward-flight 
performance of rotors with blades having as much as 3-1 taper ratio can be 
predicted with good accuracy by equations derived for uniform-chord blades, 
provided that the rotor solidity is based on the equivalent weighted 
chord Cg. 

In order to investigate the applicability of the charts to rotors 
with blades having values of twist other than -8^, the theoretical values 
of for values of twist equal to 0^, -8^, and -l6*^ were compared 

at several unstalled flight conditions. From the comparison, it appeared 
that for forward speeds ranging from the speed for minimum power to the 
maximum speeds of present-day helicopters (that is, for values of tip- 
speed ratio between approximately 0.05 and O. 50 ) the effects of twist on 
the profile power are small, particularly when considered as a percentage 
of the total power required. The importance of twist, however, is not 
primarily its effect on profile power but in the delay of stall. The 
effect of twist on stall limits is discussed in a later section of this 
paper. 


Airfoil Section Characteristics 

The three-term drag polar used in the preparation of the charts (see 
section entitled "Method of Analysis") is considered as representative of 
practical construction blades of conventional airfoil section having 
fairly accurate leading-edge profiles and rigid surfaces. The charts may 
be applied, however, to rough or poorly built blades of conventional sec- 
tion by multiplying the profile-drag — thrust ratio obtained from the 
charts by a constant "roughness" factor equal to the ratio of the average 
of the ordinates of the drag curve of the actual blade to the average of 
the ordinates of the drag curve used in the charts. If the drag curves 
do not have similar shapes, the determination of this factor should take 
into accoiKit the relative importance of different angles of attackj a 
basis for doing this by a method of "weighting" curves is discussed in 
reference 7* The. angle of attack at which stall occurs will also be 
affected by the roughness of the blade surface, and consideration should 
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be given to the surface condition when estimating the limits of validity 
of the theory. 


Stall Limits 


Satisfactory limits to the use of a theory in vhich stall is not 
considered are^ for powered flighty the conditions at which the tip of 
the retreating blade reaches its stalling angle of attack, as shown in 
references 2 and 6. For the autorotative case, limits to the theory are 
shown to consist of the conditions at which the velocity of the blade 
elements of the stalled inboard sections reach high enough values so that 
the contributions of these elements to the total thrust and torque of the 
rotor become significant. Therefore, following the procedure of previous 
MCA rotor papers (such as ref. 2), there are included on the charts of 
this paper two sets of limit lines. One set corresponds to conditions at 
which a blade element at an aximuth angle of 270^ with a relative velocity 
equal to 0.^ of the tip speed reaches angles of attack of 12^ and l6*^, 
whereas the other set corresponds to conditions at which the blade tip 
at an azimuth angle of 270^ reaches angles of attack of 12^ and l6^. 

These limit lines are designated by the symbols ay 


and a 


(1.0)(2T0°)’ respectively. 


‘(uj = 0.t)(270°) 

The 12° and l6° lines represent a range 


of angles of attack in which conventional blade airfoils would be expected 
to stall. Also, since vibration and control limitations brought on by 
blade stall occur, in general, when the calculated stall angle is exceeded 
by about 4°, the difference between the 12° and l6° lines should also be 
useful in estimating the limits to practical operating conditions of a 
rotor. Moderate amounts of stall can be approximately accounted for by 
empirical corrections to the profile power when the limit lines on the 
charts are exceeded. The basis on which these corrections may be made is 
discussed in reference 8 and the procedure is simmarized in reference k 

(pp, 266-267) • 


Theory indicates, and flight measurements have shown, that blade 
twist is effective in delaying stall. Twisting the blade so as to lower 
the pitch at the tip with respect to the pitch at the root tends to dis- 
tribute the lift more evenly along the blades and therefore minimizes the 
high angles of attack in the tip region. Combinations of Cp/Crp 
and 2CT/cra for which blade angles of attack at the specified station 
reach 12° and l6° for 0^, -8°, and -l6° twist are plotted in figure 5- 
As would be expected, these plots show that higher values of C^/a can 
be attained with negative twist before retreating blade stall is encoun- 
tered. Conversely, the greater the negative twist, the higher the tip- 
speed ratio that can be reached at a given Cfp/a before the onset of 
stall. 
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It should be noted that negative values of twist tend to decrease 
the angle of attack at the tip of the advancing blade. The advancing- 
blade-tip angle of attack is shown in figure 6 as a function of 2C^/aa 

and [I at several power conditions (as represented by the pitch values) 
for twists of -8^ and -l6^. Although the large negative angles of attack 
at the advancing -blade tip will adversely affect the performance^ this 
effect is believed to be of less importance than the benefits achieved 
by the delay in retreating blade stall. There is the possibility^ how- 
ever, that high negative advancing-blade-tip angles of attack would result 
in a contribution to blade stresses which should be considered for indi- 
vidual designs. It should be noted, however, that these tip angles were 
calculated on the basis of uniform inflow velocity, and the local upwash 
which tends to occur at the advancing tip should result in less negative 
values . 


Compressibility Limits 

The section lift and drag coefficients used in the preparation of 
the charts of this paper do not vary with Mach number. It is expected 
that the primary effect of such variation would be an increase in the 
profile-drag power if the drag-divergence Mach number were approached 
or exceeded. Therefore, the charts underestimate the power required for 
a rotor operating within the range where compressibility effects are 
encountered. It is hoped that power losses due to compressibility may 
be taken into account by adding corrections to the charts in a manner 
similar to that done for the effects of stall. The corrections probably 
could be based on results of strip analyses or on experimental data. The 
operational limits imposed by Mach number, however, are yet to be 
determined. 


CONCLUDING REMARKS 


Charts based on rotor theory have been presented from which the 
profile-drag — thrust ratio of a rotor can be determined for various com- 
binations of pitch angle, ratio of thrust coefficient to solidity, tip- 
speed ratio, and power input. The equations on which the charts are based 
have taken into account blade stall in the reversed-velocity region and 
are not limited by small-angle assumptions for blade pitch and inflow 
angles. For these reasons the method is believed to be more accurate 
than previous methods for cases wherein the rotor inflow velocity is 
relatively large, for rotors operating at steep rates of climb or descent, 
for flight at high tip-speed ratios, or for convertiplane transition 
attitudes. 
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In addition to providing a convenient means for quickly estimating 
rotor performance, the charts should he useful as a means for estimating 
the effects of changes in design variables and as a base to which correc- 
tions may be applied for the effects of stall and con5)resstbllity . Limit 
lines which indicate the stall condition of the rotor and 'vdiich serve to 
indicate the limits to practical rotor operating conditions are shown in 
the charts. 

The method of using the charts for performance estimation is outlined 
and Illustrated through computation of a sample problem. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 1, 195^. 
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Figure 1.- Concluded. 
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(a) p. = 0.05- 

Figure 2.- Thrust-coefficient— solidity ratio plotted against inflow 
ratio for blades having -8° twist. 
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Figure 2.- Concluded. 
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(a) p. = 0.05. 

Figure 3.- Profile-drag — thrust ratio for blades having -8° twist. 
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(c) ^x = 0.15. . 

Figure 5.- Continued. 
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Fig^re 4.- 


(a) Climb (7 positive). 

Cp 

Climb parameter £££_2i — 9. as a function of flight-path angle 

M- Cqi 

and parasite-drag parameter. 
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(b) Descent (7 negative). 

Figure 4 -.- Concluded. 
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